AFIT/GAE/ENY/91D-27 


ANALYSIS OF AN ADVANCED FIGHTER AIRCRAFT 
USING JET FLAP TECHNIQUES AND 
THE VORTEX LATTICE METHOD 

Thesis 

Steven P. Snyder, Captain, USAF 


AFIT/GAE/ENY/9ID-27 


92-00100 

l 



OTIC 

|£i-.ECTF. i 
>. J AN 3 8 19821 


Approved for public release; distribution unlimited 





AFIT/GAE/ENY/9ID-27 


ANALYSIS OF AN ADVANCED FIGHTER AIRCRAFT USING JET FLAP 
TECHNIQUES AND THE VORTEX LATTICE METHOD 


Thesis 

Presented to the Faculty of the School of Engineering 
of the Air Force Institute of Technology 
Air University 

In Partial Fulfillment of the 
Requirement for the Degree of 
Master of Science in Aeronautical Engineering 


Steven P. Snyder, B.S. 
Captain, USAF 

December 1991 


Approved for public release: distribution unlimited 



Acknowledgemen ts 


I would like to acknowledge Dr. Milton Franke, advisor, 
who provided support and guidance for this project: Mr. Ray 
Fredette. the designer of many interesting aircraft, who 
contributed a subject which made this study more than just 
an academic exercise: and Mr. Jim Snyder, aerodynamicist. 
who provided his personal analysis program and patiently 
endured countless questions and blank stares. I would also 
like to thank my father, Don Snyder. He sowed the seeds 
early in my life which instilled within me a love of 
aviation. Without him, I never would have made it to this 
point. 

Most of all, I would like to thank my wife. Vicke, who 
provided the endless love, support, and understanding which 
made this all bearable. 

Steven P. Snyder 



_Acoasei on F or 


NTT.; 

dt:c 










Table of Contents 

Acknowledgements . ii 

List of Figures. v 

List of Tables. xi 

List of Symbo Is.x i i 

Abstract. xv 

I . Introduction. 1 

II . Background. 4 

The Jet Flap. 4 

Applications of Blowing . 6 

Jet Flap Theory. 7 

The Vortex Lattice Method . 10 

Vortex Lattice Theory . II 

Calculation of Aerodynamic Coefficients ... 16 

Application Of The Jet Flap. 1 7 

Panel Distribution. 19 

Jet Sheet Convergence . 21 

Planform Study . 21 

Vortex Lattice Software . 22 

Data Col lect ion. 24 

III. Aircraft Configuration . 25 

Design Philosophy . 25 

Requirements and Figures of Merit . 26 

The Computational Model . 28 

Pane 1 1 ing. 28 

Blowing Coefficient . JO 

Nozzle Vector Angle . 31 

Modification for Control Surface Sizing ... 31 

IV. Verification. 3 2 

Verification of Vortex Lattice Code . 32 

Two Dimensional Case With Camber (NACA 

64-206). 32 

Jet Sheet Convergence . 32 

Three Dimensional Wing With Camber. 33 

NADC Jet Flapped Fighter Study. 33 

V. Resu Its. 38 

Force Accounting. 38 

b.10JF-B ResuIts . 40 

2.1 Influence of Blowing Coefficient .... 40 

2.2 Influence of Flap Deflection. 44 

i i i 







































Influence of Nozzle Width . 49 

Stability and Control Results . 51 

Planform'Study Results . 52 

VI. Conclusions and Recommendations . 56 

Bibliography . 58 

Appendix A: Philosophy of Designer . 60 

Appendix B: b.lOJF-B Aerodynamic Data . 63 

Appendix C: Stability and Control Data. 86 

Appendix D: Planform Data. 92 

Vita. 122 


i v 
















Li st of Figures 


Figure 1 Wright Lab b.lOJF-B Two View. 2 

Figure 2 Jet Control Volume (Adapted From Spence) . 8 

Figure 3 Vortex Filament (From Bertin. et al.) ... 13 

Figure 4 Horseshoe Vortex (From Bertin, et al.) ... 14 

Figure 5 Paneled Wing. 14 

Figure 6 b.lOJF-B Panel Model . 28 

Figure 7 NASA Transonic, Sharp Leading Edged Airfoil . 29 

Figure 8 DTRC Fighter Comparison Lift Curves . 36 

Figure 9 DTRC Fighter Comparison Drag Data. 36 

Figure 10 DTRC Fighter Comparison Moment Data .... 37 

Figure 11 C Djet Schemat ic. 40 

Figure 12 Variation in C La with Respect to Cj. 42 

Figure 13 Example b.lOJF-B Lift Data, 8j=40°. 42 

Figure 14 Variation of C L8j with Respect to Cj .... 43 

Figure 15 Example b.lOJF-B Moment Data, 83 = 40 " .... 44 

Figure 16 Variation in Static Margin with Respect to 

Cj. 45 

Figure 17 Example b.lOJF-B Drag Data, 8j=40’ 46 

Figure 18 Example b.lOJF-B C L vs 83 , BC=2 . 47 

Figure 19 C La Var iat ion with Respect to 8 j. 4 "? 

Figure 20 Example b.lOJF-B Drag Data, BC=2 . 48 

Figure 21 Example b.lOJF-B Moment Data. BC=2 . 49 

Figure 22 Variation in Static Margin with Respect to 

6 j. 49 

v 





















Figure 23 Example b.lOJF-B Nozzle Width Comparison, 

Lift . 50 

Figure 24 Example b.lOJF-B Nozzle Width Comparison, 

Drag. 51 

Figure 25 Example b.lOJF-B Nozzle Width Comparison, 

Moment. 5 1 

Figure 26 Example b.lOJF-B With Canard Removed, 

Moment. 53 

Figure 27 Lift Curve Variation with Respect to AR . . 54 

Figure 28 Drag Polar Variation with Respect to AR . . 54 

Figure 29 Moment Variation with Respect to AR .... 55 

Figure 30 Variation of Lift with Respect to A, ... . 55 

Figure 31 Variation in Drag with Respect to A ... . 56 

Figure 32 Moment Variation with Respect to A .... 56 

Figure 33 Lift Curve for 8j=-10° 65 

Figure 34 Drag Polar for Sj=-10° 65 

Figure 35 Pitching Moment About the CG for 6 j=-10° . . 66 

Figure 36 Lift Variation Due to Cj for Sj=-10°. G 6 

Figure 37 Lift Curve for Sj=0° 67 

Figure 38 Drag Polar for 8j=0° 67 

Figure 39 Pitching Moment About the CG for 8j=0° ... 68 

Figure 40 Lift Variation Due to Cj for 8j=0°. 68 

Figure 41 Lift Curve for 8-10°. 69 

Figure 42 Drag Polar for 83 = 10 °. 69 

Figure 43 Pitching Moment About the CG for 8=10° ... 70 

vi 


















Figure 44 Lift Variation Due to Cj for 6j=10 o . 70 

Figure 45 Lift Curve for 6j= 20°. 71 

Figure 46 Drag Polar for 6j=20°. 71 

Figure 47 Pitching Moment About the CG for 8j=20° ... 72 

Figure 48 Lift Variation Due to Cj for 8j=20°. 72 

Figure 49 Lift Curve for 8^30°. 7 3 

Figure 50 Drag Polar for 6j=30°. 73 

Figure 51 Pitching Moment About the CG for 6j=30“ ... 74 

Figure 52 Lift Variation Due to Cj for 8j=30°. 74 

Figure 53 Lift Curve for 8^=40'. 75 

Figure 54 Drag Polar for 6^=40°. 75 

Figure 55 Pitching Moment About the CG for 5j=40° ... 76 

Figure 56 Lift Variation Due to Cj for 6j=40°. 76 

Figure 57 Lift Curve for BC=0. 7 7 

Figure 58 Drag Polar for BC=0. 77 

Figure 59 Pitching Moment About the CG for BC=0 .... 78 

Figure 60 Lift Variation Due to 8j for BC=0. 78 

Figure 61 Lift Curve for BC=1. 79 

Figure 62 Drag Polar for BC=1. 7 9 

Figure 63 Pitching Moment About the CG for BC=1 .... 80 

Figure 64 Lift Variation Due to 8j for BC=1. 80 

Figure 65 Lift Curve for BC = 2. 81 

Figure 66 Drag Polar for BC = 2. 81 

Figure 67 Pitching Moment About the CG for BC=2 .... 82 

Figure 68 Lift Variation Due to 6j for BC=2. 82 


v i i 

























Figure 69 Lift Curve for BC = 3. 83 

Figure 70 Drag Polar for BC=3. 83 

Figure 71 Pitching Moment About the CG for BC=3 .... 84 

Figure 72 Lift Variation Due to Sj for BC=3. 84 

Figure 73 Lift Curve for BC=4. 85 

Figure 74 Drag Polar for BC=4. 85 

Figure 75 Pitching Moment About the CG for BC=4 .... 86 

Figure 76 Lift Variation Due to 8j for BC=4. 86 

Figure 77 Pitching Moment About CG, 6j=-10°, No Canard 88 

Figure 78 Pitching Moment About CG, Sj=0*, No Canard . 88 

Figure 79 Pitching Moment About CG, 8j=10*, No Canard . 89 

Figure 80 Pitching Moment About CG, 6j=20°, No Canard . 89 

Figure 81 Pitching Moment About CG, 8j=30°, No Canard . 90 

Figure 82 Pitching Moment About CG, Sj=40°, No Canard . 90 

Figure 83 Pitching Moment About CG, BC=0, No Canard . . 91 

Figure 84 Pitching Moment About CG, BC=1, No Canard . . 91 

Figure 85 Pitching Moment About CG, BC=2, No Canard . . 92 

Figure 86 Lift Curve Comparison for Sj, AR=2, >1=0.05 . 94 

Figure 87 Lift Curve Comparison for BC, AR=2, A.=0.05 . 94 

Figure 88 Drag Polar Comparison for 8j, AR=2, A,=0.05 . 95 

Figure 89 Drag Polar Comparison for BC, AR=2, A. = 0.05 . 95 

Figure 90 Moment Comparison for 8j, AR=2, >1=0.05 ... 96 

Figure 91 Moment Comparison for BC, AR=2, A.=0.05 ... 96 

Figure 92 Lift vs 8j, AR=2, A, = 0.05 . 9" 

Figure 93 Lift vs Cj, AR=2, A,=0.05 . 9' 


vi i i 














Figure 94 Lift Curve Comparison for Sj, AR=2, >=0.25 . 98 

Figure 95 Lift' Curve Comparison for BC, AR=2, >=0.25 . 98 

Figure 96 Drag Polar Comparison for 8j, AR = 2, >=0.25 . 99 

Figure 97 Drag Polar Comparison for BC, AR=2, >=0.25 . 99 

Figure 98 Moment Comparison for 8j, AR=2, >=0.25 . . . 100 

Figure 99 Moment Comparison for BC, AR=2, >1 = 0.25 . . . 100 

Figure 100 Lift vs Sj, AR=2, >1=0.25. 101 

Figure 101 Lift vs Cj, AR=2, >=0.25. 101 

Figure 102 Lift Curve Comparison for 8j, AR=3, >=0.05 . 102 

Figure 103 Lift Curve Comparison for BC, AR=3, >1 = 0.05 . 102 

Figure 104 Drag Polar Comparison for 8j, AR=3, >=0.05 . 103 

Figure 105 Drag Polar Comparison for BC, AR=3, >=0.05 . 103 

Figure 106 Moment Comparison for 8j, AR=3, >=0.05 . . . 104 

Figure 107 Moment Comparison for BC, AR=3, >.=0.05 . . . 104 

Figure 108 Lift vs 8j, AR=3, >,=0.05. 105 

Figure 109 Lift vs Cj, AR=3, >. = 0.05. 105 

Figure 110 Lift Curve Comparison for 6j, AR=3, >1 = 0.25 . 106 

Figure 111 Lift Curve Comparison for BC, AR=3, >1=0.25 . 106 

Figure 112 Drag Polar Comparison for 8j, AR=3, >1 = 0.25 . 107 

Figure 113 Drag Polar Comparison for BC, AR=3, >,=0.25 . 107 

Figure 114 Moment Comparison for 8j, AR=3, >=0.25 . . . 108 

Figure 115 Moment Comparison for BC, AR=3, >=0.25 . . . 108 

Figure 116 Lift vs 8 jt AR=3, >=0.25. 109 

Figure 117 Lift vs Cj, AR=3, > = 0.25. 109 

Figure 118 Lift Curve Comparison for 6j, AR=3, >=0.5 . 110 

ix 













Figure 119 Lift Curve Comparison for BC, AR = 3, A, = 0.5 . 110 

Figure 120 Drag Polar Comparison for Sj, AR=3, X = 0.5 . Ill 

Figure 121 Drag Polar Comparison for BC, AR=3, A,=0.5 . Ill 

Figure 122 Moment Comparison for 6j, AR=3, >1=0.5 . . . 112 

Figure 123 Moment Comparison for BC, AR=3, A.=0.5 . . . 112 

Figure 124 Lift vs 8 j5 AR=3, A = 0.5 113 

Figure 125 Lift vs Cj, AR=3, X = 0.5 113 

Figure 126 Lift Curve Comparison for 8j, AR=5, >1=0.5 . 114 

Figure 127 Lift Curve Comparison for BC, AR=5, A=0.5 . 114 

Figure 128 Drag Polar Comparison for 8j, AR=5, X=0.5 . 115 

Figure 129 Drag Polar Comparison for BC, AR=5, >1=0.5 . 115 

Figure 130 Moment Comparison for Sj, AR=5, X=0.5 . . . 116 

Figure 131 Moment Comparison for BC, AR=5, >1=0.5 . . . 116 

Figure 132 Lift vs Sj, AR=5, X=0.5 117 

Figure 133 Lift vs Cj, AR=5, X = 0. 5 117 

Figure 134 Lift Curve Comparison for 8j, AR=5, il=l . . 118 

Figure 135 Lift Curve Comparison for BC, AR=5, >1=1 . . 118 

Figure 136 Drag Polar Comparison for 6j, AR=5, X=1 . . 119 

Figure 137 Drag Polar Comparison for BC, AR=5, >1=1 . . 119 

Figure 138 Moment Comparison for 8j, AR=5, A=1 .... 120 

Figure 139 Moment Comparison for BC, AR=5, A=1 .... 120 

Figure 140 Lift vs 8j, AR=5 , >1=1 121 

Figure 141 Lift vs Cj, AR=5, A=1 121 

x 













List of Tables 


Table I Planform Study Wing Parameters. 22 

Table II b.lOJF-B Engine Characteristics . 30 

Table III LATJET Comparison for N'ACA 64-206 32 

Table IV Analytic and Computational Values for 

b. 10JF-B Wing. 33 

Table V b.lOJF-B Constant Parameters . 41 

Table VI b.lOJF-B Blowing Conditions . 42 

Table VII b.lOJF-B Dimensions and Weight. 62 













List of Symbols 


AR 

b 

BC 

c 

c 

Coi 

eg 

Cj 

C L a 

c, 

Cl 

c„ 

C.,„ 

C, 

deg 

F’ 

FS 

1 

L 

L’ 

m j 

m’ 

N 

P 


Aspect Ratio 
Span 

Blowing Condition 
Chord 

Mean Aerodynamic Chord 
Induced Drag Coefficient 
Center of Gravity 
Blowing Coefficient, Gross 
Lift Curve Slope 

Two Dimensional Lift Coefficient 

Three Dimensional Lift Coefficient 

Moment Coefficient 

Influence Coefficient 

Blowing Coefficient, 2D 

Degree 

Force Per Unit Span Acting on Jet 

Fuselage Station 

length 

Lift 

Lift Per Unit Span 

Mass Flow Rate of Jet 

Mass Flow Rate of Jet Per Unit Span 

Total Number of Panels 

Pressure 

x i i 






q 

r 

R 

S 

SM 

TR 

V 
v j 

a 

Y 

r c 

r j 

r -n 

8.n 

0 

X 

A c /2 

p 

Subscript s 

CLmax 


j 

L = 0 


Dynamic Pressure 

Radial Distance 

Jet Sheet Deflection Radius 

Area 

Static Margin 
Taper Ratio 
Induced Velocity 
Jet Velocity 

Nondimentionalized Downwash 
Angle of attack 
Circulation per Unit Area 
Wing Circulation Per Unit Span 
Jet Circulation Per Unit Span 
Local Circulation 
Flap Deflection Angle 
Kronecker Delta 
Jet Sheet Curvature Angle 
Taper Ratio 

Sweep of Wing Half Chord Line 
Air Density 

Pertaining Maximum Lift 

Induced Parameter 

Pertaining to the Jet 

Pertaining to Zero Lift Conditions 

x i i i 







f Current Interest 


icing Panel 
earn Condition 






Abstract 


A computer code employing the vortex lattice method 
with modified boundary conditions was used to determine the 
induced aerodynamic characteristics of high aspect ratio, 
vectoring, exhaust nozzles located at the wing root of a 
canard configured fighter aircraft. Comparison with 
existing wind tunnel data verified results of the method. 

The exhaust was modeled as a singularly blown jet flap at 
deflection angles of -10, 0, 10, 20, 30, and 40 degrees. 

Jet momentum coefficients were based gross engine thrust for 
maximum afterburner and military power settings at a Mach 
number of 0.6 and an altitude of 20,000 ft. Lift, induced 
drag, and pitching moment coefficients were calculated for 
untrimmed conditions. The nozzle provided lift augmentation 
at all deflections and blowing conditions, a reduction in 
induced drag at high lift coefficients, and an increased 
nose down pitching moment. An optimum flap deflection to 
achieve minimum induced drag existed for each unique lift 
and blowing condition. Measurement of the static margin 
showed that the aircraft as configured was statically 
unstable. Example cases of comparable wing planforms 
varying in aspect ratio from 2 to 5 and taper ratio from 
0.05 to 1 showed similar behavior. 
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ANALYSIS OF AN ADVANCED FIGHTER AIRCRAFT USING JET FLAP 


TECHNIQUES AND THE VORTEX LATTICE METHOD 


I. Introduction 


Conceptual aircraft design generally concerns itself 
with trends or technology that may be 15 to 20 years in the 
future. At this early point in the life of an aircraft, the 
designer must be able to comprehend how proposed 
technologies will enable accomplishment of requirements put 
forth. This comprehension leads to well founded decisions 
throughout the process. One technology which has been 
proposed for use on advanced aircraft for many years is the 
jet flap. 

In the case of this study, the designer, Mr Raymond 
Fredette, proposed the use of a very high aspect ratio jet 
exhaust nozzle located at the wing/body juncture of a canard 
fighter configurat ion.(See Figure 1) One reason for 
incorporating this type of nozzle was to improve the drag 
polar of the aircraft due to the increase in circulation 
provided by a jet flap. There were also survivability 
benefits and weight penalties to be investigated. The key 
question was whether or not the benefits would outweigh the 


pena 11 ies. 















Fredette had no verified aerodynamic methods to predict 
the effects of the nozzle. Additionally, the aircraft 
itself was pushing the capability limits of existing 
empirical methods because of its unusual configuration. He 
believed that some sort of computational method was needed 
to generate useful aerodynamic data. By treating the thin, 
vectored exhaust like a jet flap, it was felt that 
worthwhile aerodynamic coefficients could be determined. 
After studying several methods of calculating jet flap 
effects, a vortex lattice method with modified boundary 
conditions was selected. 

A secondary study of generic wing planforms similar in 
aspect ratio and taper ratio to the b.lOJF-B wing was also 
performed in order to demonstrate the effects of varying 
wing parameters on the aerodynamic characteristics of a jet 
flapped wing. 
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I I. Background 


1 The Jet Flap 

According to Kohlman, (Ref 2), Shubauer first proposed 
the jet flap in 1932. Kohlman goes on to state that 
Hagedorn and Ruden carried out the initial experimental work 
in 1938. The concept involves a sheet of high velocity air 
blown from a thin slot near the trailing edge of a wing. 
Because the sheet has a finite amount of momentum, it can 
support a pressure difference which deflects the freestream. 
(Ref 14) This is analogous to the Kutta condition being met 
somewhere aft of the trailing edge. The result is an 
increase in lift beyond the circulation lift of the wing and 
the lift component of jet thrust. According to McCormick, 
under strong blowing conditions, the jet can also act like a 
boundary layer control device and delay leading edge 
separation, resulting in 6 to 8 degree increments in a CL-8Jc . 

The change in induced drag compared to a standard wing 
depends on the amount of blowing and flap deflection, 6j. 
Induced drag reduction is realized at high lift conditions 
while some penalty is paid at low lift. There is a point in 
between these conditions where the unblown and blown drag 
polars cross over each other. This point is dependant on 
the amount of blowing and the flap deflection. 

Theoretically, all thrust from the deflected jet is 
recovered if it is turned parallel to the freestream. 
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However, due to a finite value of downwash at infinity, a 
fraction remains unrecovered and is accounted for as an 
induced drag component, C Djet . Some thrust is also lost due 
to friction at the jet/freestream boundary.(Ref 14) One 
penalty paid is in the form of a severe nose down pitching 
moment at high lift conditions. The drag associated with 
trimming this moment can be significant enough to outweigh 
many benefits. 

Jet flaps come in two distinct varieties. The first is 
a singularly blown flap where the jet sheet exits the wing 
at the trailing edge and is vectored by a nozzle 
arrangement. The second type is known as a blown flap. It 
blows the jet over a mechanical flap which deflects it due 
to the Coanda effect. Each concept has its shortcomings. 

The singularly blown flap is tough to incorporate both 
structurally and mechanically into a thin trailing edge, 
while the blown flap has problems with the hot exhaust 
interacting with the flap. This study will model the jet 
exhaust as a singularly blown flap. 

Confusion sometimes exists between the concepts of the 
jet flap and circulation control since they both use a form 
of blowing to achieve high lift conditions. With 
circulation control, a thin jet of air is blown over a 
rounded trailing edge, eliminating the Kutta condition. 

This moves the aft stagnation point to the lower surface. 
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thereby increasing circulation and lift.(Ref 14) On the 
other hand, the primary purpose of blowing in a jet flap is 
to add momentum to the wake in order to support a pressure 
discontinuity. Each method has different means of achieving 
similar objectives. 

1.1 Applications of Blowing 

One of the few purely jet flapped aircraft to fly, was 
the Hunting H.126. (Ref 2, 15) Engine air was ducted solely 
to full span, blown flaps along the trailing edge of this 
research aircraft. All thrust used to propel the aircraft 
in forward flight came from the jet flaps. 

Two forms of blown flaps were utilized in competing 
prototypes for the Air Force Advanced Medium STOL Transport 
(AMST) program in the 1970s. (Ref 15) McDonnell Douglas 
used an externally blown flap on its entry, the YC-15. The 
exhaust of its underwing mounted turbofan engines was blown 
through large two-segment slotted flaps. Boeing employed 
upper surface blowing on the YC-14. The exhaust from its 
over the wing mounted engines was deflected downward over 
large flaps via the Coanda effect. 

The Navy tested a rather bulky circulation control 
high lift device on a Grumman A-6 in the late ’70s. (Ref 2) 
Bleed air pulled from the engines was ducted through plenums 
mounted at the trailing edges. A slot running lengthwise 
along the plenums provided an exit for the jet sheet. The 
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rounded outer surface of the plenum acted as a Coanda 
surface in order to deflect the flow. In order to counter 
the strong nose down pitching moment, the horizontal 
stabilizer was increased in area and recountoured with an 
extreme negative camber. 

Blowing is typically associated with lift enhancement 
in aircraft. However, reference 6 proposed the use of a 
circulation control stern plane on a submarine. The 
purpose was to prevent a catastrophic crash dive which could 
be caused by a mechanical stern plane becoming jammed. 

1.2 Jet Flap Theory 

The basis for most analytical theory of the jet flap 
can be attributed to D. A. Spence of the Royal Aircraft 
Establishment. (Ref 17) The following development was 
adapted from Reference 14. Spence used a control volume 
approach to represent the jet exiting the trailing edge of 
the wing. Although it entered at some angle against the 
free stream, he reasoned that it would eventually be turned 
parallel. The pressure, Ap which turned the jet. 
transferred itself in part to the lifting surface. Figure 2 
depicts this concept in schematic form. 






r 

i 

i 

i 

z* 


Figure 2 

Applying the momentum equation to the control volume, 
m'jVj A0 = ApJ?A0 ( 1) 

and solving for Ap. 

( 2 ) 

The force per unit width acting on the jet is, 

F* = Api?A0 (3) 

Using a vortex with strength per unit width, the force can 
be represented by: 

F* = p m VjfRA9 (4) 

Solving for y as circulation due to the jet, yj, 
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J « + *j 

< w, 


I p+Ap 

Jet Control Volume (Adapted From Spence) 



With a large radius of curvature and small A0 , this becomes 


yj = 


-m'jVj d 2 z 

P.V m dx 2 


( 6 ) 


Integrating Yj from the trailing edge to infinity yields Tj, 
jet induced circulation per unit span. 


* 


' m j v j r cPz 
P.v. J dx 2 


dx « 


(7) 


given small values of a and 6j. The total two dimentional 
lift coefficient is then expressed as a combination of 
factors due to the airfoil and the jet. 


2T C iriiV , 

~V~c 


qjc 3 


2T C 


+ C («+*J 


(8) 


where the blowing coefficient, C„, is the primary 
nondimensiona1 parameter used in describing jet momentum in 
two dimentions. Cj is a gross blowing, or jet coefficient 
for three dimentions. 


<L,S 


(9) 


The solution of T c and fj involves a pair of 
simultaneously solved integro-differentia 1 equations. 






Spence graciously provided curve fit solutions to the 
equations, Eq ( 10 ) . 


acj 
~d6 = 


4.:C ll + 0.151C„ 2 +0.139 CJ 


dC 2 

da 


2ltll+0.151C* +0.219 CJ 


„ dc, dc 2 - 

1 a« m 


( 10 ) 


These allow painless determination of jet flapped airfoil 
solutions. 

Spence continued on to a three dimensional solution 
with Maslcel. A very good description of this method is 
found in Reference 14. However, the solution assumes that 
the lift distribution and the blowing coefficient vary 
elliptically over the span. Because this is an idealized 
case, and the objective of this study was to analyze a 
specific configuration, a method with more flexibility was 
pursued. 


2 The Vortex Lattice Method 

One of the primary computational methods used to model 
the jet flap has been the vortex lattice method. Using 
boundary conditions similar to those set up by Spence, the 
jet sheet can be modeled as a component interacting with 
wing and body.(Ref 4) 
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The level of detail involved with methods varies 


greatly from code to code. The author wrote a relatively 
short BASIC program as an undergraduate to solve for the 
incompressible aerodynamic coefficients of basic planar 
geometries. On the other hand, Douglas Aircraft developed 
the Elementary Vortex Distribution, or EVD, panelling code 
which modeled circulation singularities at the wing leading 
edge with square root singularities and hinge lines with 
logarithmic singu1arities.(Ref 13) 

2.1 Vortex Lattice Theory 

The theoretical development described here is actually 
from a combination of sources. (Ref 3,4,5) Some variable 
and indices have been changed to assure uniformity. 
Essentially, the vortex lattice method models the velocity 
induced by a wing of finite aspect ratio with a distribution 
of horseshoe vortices. Each vortex is of a strength which 
satisfies a set of boundary conditions. These strengths are 
solved simultaneously ever the planform with a system of 
linear equations. Given the strengths, the aerodynamic 
characteristics of the configuration can then be determined. 

The foundation of vortex lattice theory is the Biot 
Savart law, which is used in determining the velocity 
induced by any given vortex filament at a particular control 
point. In vector form: 





r a (c?lxf) 


(in 


dV = 


integrating Eq (11) determines the induced velocity at a 
given point in space for a given vortex strength, r n . 


7—2- fainddd = —^2-(cos0j-cose 2 ) (12) 


Pi x f 2 | 


COS0! = 


COS0, = —-— 


(13) 


r n A x 

411 Pi X / 2 | 2 


y - * • t2 1 f ■ (h-h 

- - -.*r° ( r i 


(14) 


The most general form of Eq (14) is the formula for the 
induced velocity of a straight infinite vortex filament, 


r a 

2 nr 


(15) 


Figure 3 shows a generic vortex filament and its associated 
parameters. 

According to the Helmholtz theorem, any vortex in 
inviscid flow must either, extend to infinity, terminate at 
a wall, or form a closed loop. (Ref 14) In the case of the 
horseshoe vortex in Figure 4, points A and B define the 
bound portion of the vortex, and both ends trail off to 





positive infinity. The induced velocity 


determined once the influences of the individual segment 
are known. Due to the linear nature of potential flow ai 
its elementary flows, the influences of the legs may thei 
added together to get the total velocity. In turn, the 
induced velocity at any given point may be determined by 
summing the contributions of all the other horseshoe 
vortices. 

Assuming the wing in Figure 5 is to be studied, 
modeling begins with dividing the planform up into spanw 
and chordwise panels. Each panel corresponds to a horse 
vortex.(The distribution of these panels will be discuss 
later.) The sides of each panel are parallel to the X a 


iw the taper o 








panel and the free ends are allowed to trail back to 
infinity in the X direction. Surface boundary conditions 
are enforced at the control points located at the 3/4 chord 
location, midway between the sides. 

The following notation is helpful when summing up 
velocities. The panel number of the current control point 
is represented by the m index, and n represents a horseshoe 
vortex influencing it, such that: 

06 ) 

C is the influence coefficient based on the location and 
orientation of the vortex. It is the only known portion of 
Eq (14) since it is determined by geometry alone. Summing 
the individual velocities together, 

07) 

n«l 

Where the right hand side is the normalized velocity at 
induced at control point m. The left hand side sets up a 
system of linear equations which may be solved for r n /Vv 
Boundary conditions specify the normalized velocity. 

2.1.1 Boundary Conditions 
Boundary conditions for the vortex lattice method 
result from the requirement that the flow be tangential to 
the surface at the control points located at the 3/4 panel 
chord, and mid panel span. The easiest way to enforce this 
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constraint is to force the normal component of velocity to 
zero at all control points. 

(V m +V m )‘A = 0 (18) 

The system of linear equations is then solved using the 
matrix inversion technique of choice. 

2.2 Calculation of Aerodynamic Coefficients 
Once the individual panel circulations are calculated, 
the values may be used to determine the forces acting on 
individual panels. Starting with the Kutta-Jukowski law in 
its most general form, 

s 1 . l 

F - jT p^xTcte (19) 

a 1 

The lift per unit span becomes 

L’ - p_V_r ( 20 ) 

and C L is solved in terms of the local singularities in both 
the chordwise and spanwise directions. 



N is the total number of panels on one half of a symmetric 
wing. The moment is simply the summation of the various 
panel lift contributions multiplied by the moment arm from 
the point of interest to the panel vortex location. 
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The induced drag coefficient is slightly more complicated. 
For the near field drag value used in this study, the 
induced velocity normal to the surface at the local bound 
vortex must be determined. This velocity then takes the 
place of V a in Eq (20), and the summation becomes 


c », * jEE c «. 


r r 

— Ay„ 
1 v m v m ° 


2.3 Application Of The Jet Flap 

Applying the jet flap to the vortex lattice method is 
an exercise in enforcing a specific set of boundary 
conditions over an additional set of panels representing the 
jet sheet. Since the jet sheet supports a pressure 
discontinuity, it can be modeled with a typical horseshoe 
vortex panel. The following development was adapted from 
Bonner et al. Some variables and indices have been renamed 
for the sake of continuity. 

The change in jet inclination due to the pressure 
discontinuity causes the major difference to boundary 
conditions. 

First, starting over with Eq (4) 



F* = p m VjfRAd 


( 4 ) 


and note that for a large radius of curvature, 

A is d 2 z _ dw 
R dx 2 ~ dx 


(24) 


where w represents the downwash nondimentionaIized with 
respect to V a and dw/dx represents its change along the 
length of the jet. Then 


ll = m 'l v i = = 1 c dw 

v - p m vtR p m Vt dx 2 * dx 


(25) 


or 




(26) 


Integrating this equation in the x direction, between 
control points 


cC - / -f '- 2 / • 0 


(27) 


resulting in 


cC^fcV'^i] 


V m 


(28) 


Given the fact that w is a function of T, the governing 
equation becomes 
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( 29 ) 


n»l v m 

where 6 -n is the Kronecker delta. (Not to be confused with 
6j, the flap deflection angle) The jet panel adjacent to 
the jet exit must take into account the flap deflection 
angle relative to the wing. 

cC it [w m -(W„. 1 +b j ))-2^ = 0 (30) 

Now. the summation becomes: 

X) [ c ^vfPm.n~ C m-i.n) _2 ®«,n] ~y = C( ~v^i (31) 

n-l r « 

Essentially, this means that the influence coefficient on 
the jet, C.n, becomes a function of the traditional 
influence coefficient on the current panel and the one 
immediately aft of it. 

- CC„(C.. 0 -C M .J-2S„.„ 

on the jet, and the right hand side becomes: 

cCfij on the jet, adjecent to the nozzle 
0 elsewhere 

2.4 Panel Distribution 
Panel distribution remains as one of the few black arts 
left regarding paneling methods. The fact that the vortex 


(32) 


(33) 
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lattice method actually solves the flow on 1v at the control 
points implies that it probably misrepresents certain 
phenomena. One of the primary things it does miss are the 
singularities in the vorticity distribution at the leading 
edge and hinge lines. One solution to this problem 
increases panel density in the location of a singularity. 
However, this quickly eats up computer time and storage, and 
some codes may only have a limited number of panels to use. 
The McDonnell Douglas EVD Method (Ref 13) models the 
chordwise circulation at these points with logarithmic and 
inverse square root singularities. Another approach to 
increase convergence, using significantly fewer panels, is 
to space the wing tips and hingelines fractions of panel 
widths from their actual boundaries. 

The number of spanwise panels necessary for a given 
convergence has been shown to be significantly reduced by 
insetting the wing tip 1/4 of a panel span. Hough proved 
this by deduction for an elliptic planform. By induction, 
he shows it to be the case for wings with sweep and constant 
chord. 

A similar method exists when modeling the circulation 
singularity existing at a flap hingeline. The bound vortex 
of the flap panel will better represent the circulation if 
it overlaps the flap and wing by approximately 1/4 of a 
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panel chord. (Ref 19) This places it closer to the 
hinge line. 

In general, adjacent panels should be kept relatively 
similar in size and shape. Significant differences will 
exist between the root and tip panels of a highly tapered 
wing, but the variation is minimal from panel to panel. 

This rule becomes important when joining separate components 
like a wing and flap. Adjacent panels should be aligned and 
of the same span. The chords should also be similar. 

2.5 Jet Sheet Convergence 

Theoretically, the jet sheet extends to infinity. To 
limit the size of the model, an investigation was conducted 
to determine the minimum sheet length required to achieve 
acceptable convergence. A wing planform was run with jet 
sheets varying in length from 1 to 10 chord lengths. The 
lift coefficient converged to within 0.1% of the final value 
in three chord lengths. The bound vortex strengths on the 
jet sheet were negligible aft of this point. 

3 Planform Study 

This additional investigation was included to give the 
reader some idea of how the characteristics of a wing 
planform affect the influence of the jet flap, especially 
for lower aspect ratio wings. The study varied planform AR 
from 2 to 5 and X from 0.05 to 1. Table I lists the 
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parameters which specifically describe each planform. Each 
case was run for the following values of a.Oj and C^. 
a=0, 5 , 10, 15, 20 deg 

6j = 0, 15, 30 deg 

= 0, 0.5, 1.0 


Table I Planform Study Wing Parameters 


AR 

k 

b ( 1 ) 

, 

S ( 1") 

c ( 1 ) i 

A c/2 (deg) 

2 

0.05 

2.121 

2.25 

1.35 | 

0 

2 

0.25 

2. 121 

2.25 

1 

1.188 | 

0 

3 

0.05 

2.598 

2.25 

1.102 | 

0 

3 

0.25 

2.598 

2.25 

0.97 | 

0 

3 

0.5 

2.598 

2.25 

0.898 

0 

5 

0.5 

3.354 

2.25 

0.696 j 

0 

5 

1 

1 

3.354 

2.25 | 

0.671 | 

0 


4 Vortex Lattice Software 

The computer code used in this study, LATJET, was 
written by Mr James Snyder, ASD/XRH. (Ref 18) LATJET allows 
relatively complex shapes to be developed. Components are 
developed by defining chords from root to tip. 

Subcomponents between chords are divided into evenly spaced 
spanwise and chordwise panels. Up to 300 panels may be 


def ined. 











Camber and twist are modeled by modifying the local 
angle of attack at a given pane 1 to account for the 
variation in surface slopes. The local angle of attack is 
factored into the boundary conditions when the vortex 
lattice method constructs the set of linear equations. 

Slopes can be input by hand or a table may be created using 
CAMBER. CAMBER is a FORTRAN code developed in house by 
ASD/XRH for producing specific camber distributions. Camber 
lines may be defined in one of several ways.(Ref 7) One 
method generates camber lines given a combination of one or 
more NACA airfoil types. Lines are defined at both inboard 
and outboard chords so that spanwise variations can be 
determined by interpolation. 

The jet sheet is defined as an individual component 
linked to the trailing edge of another. C w and 8j are 
specified on a strip by strip basis. 

At subsonic speeds, the program uses vortex lattice 
panels, while above Mach 1, constant pressure panels are 
employed. A Prandt1-Glauret correction accounts for 
compressibility. The set of linear equations is solved 
using Lower Upper decomposition. This solution method 
allows a significant reduction in computation time if 
multiple boundary conditions such as camber, control 
deflections or a, are to be run for the same geometry. 
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The LATJET output file can vary greatly in size 
depending the settings of the user defined print options. 
Individual component, strip or panel information may be 
output in addition to a final table which contains the 
resulting aerodynamic, moment and jet flap coefficients of 
the model being studied. 

Since the method is predicated upon potential flow, no 
viscous effects are modeled. This is especially important 
when looking at the C L vs a curves. No maximum lift 
coefficient is implied. The relationships are linear. In 
this case, the designer will use empirical methods to 
determine the effects of viscosity. All components are 
assumed to be thin, t/c << 1. 

5 Data Collect ion 

LATJET was run on the VMS cluster at the Air Force 
Institute of Technology. Individual data cases were piggy¬ 
backed into large files with similar blowing conditions. 
These files were then processed in the background mode using 
multiple batch files. A FORTRAN sorting code distilled the 
large output files into manageable sized ASCII files. 
Finally, these data files were imported into Quattro® Pro 
3.01 where they were parsed and lumped into a large 
relational data-base. Data could be plotted at will for any 
available choice of parameters using a custom made 


spreadsheet. 





Aircraft Configuration 


111 


1 Design Philosophy 

The personal philosophy behind the design of the 
b.lOJF-B is described by Fredette in Appendix A and 
reference 8. One potential benefit derived from choosing 
the high aspect ratio nozzle is survivability. It is widely 
circulated in open literature that a nozzle of this type 
masks hot parts and mixes the exhaust plume.(Ref 16) In 
addition to the jet flap benefits mentioned, there is also 
the prospect that this type of nozzle can reduce the boat- 
tail drag (drag caused by the separation of airflow over the 
aft portion of the fuselage), if it can induce favorable 
flow conditions over the aft portion of the fuselage. 

The advantages of the jet flap seem very lucrative. 
However, there are several disadvantages. The primary are: 
nozzle weight and the additional fuselage weight acquired 
when accommodating the greater transition section weight, 
length, and width of the nozzle. Using first order 
estimates, the nozzl% and transition section combination can 
constitute from 40 to 60 percent of the bare engine 
weight.(Ref 8) Weight penalties magnify several times when 
actual take off gross weight is addressed. In the case of 
the b.lOJF-B configuration, each additional pound of nozzle 
weight increased the gross weight by 3.9 pounds when 
considering the extra body weight and additional fuel 
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required to meet mission radius requirements. Complexity is 
another detractor. However, standard nozzles, especially 
two dimensional vectoring types, are approaching the jet 
flap in overall complexity. The extreme length of the 
nozzle and the fact that it must transition the flow from a 
round to rectangular cross section indicates that propulsive 
losses will be higher than standard. If thv. nozzle reduces 
boat-tail drag, the length could be shortened. The shorter 
transition will create steeper slopes at the aft end, but 
separation conditions might be avoided due to favorable 
influence tlve flowfield by the jet. 

2 Requirements and Figures of Merit 

When determining whether or not the nozzle will be of 
benefit to the aircraft, there are several performance 
requirements to be considered. These requirements, set at 
the beginning of the design study, can be divided into two 
main categories, maneuverability and mission performance. 
Maneuverability concerns itself with the aircraft’s 
performance at specific points in the operating envelope 
usually defined with Mach number, altitude, and power 
setting. The different criteria include instantaneous turn 
rate, sustained turn rate and specific excess power. 

Mission performance deals with ranges, loiter times, take¬ 
off d i stance, etc. 
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Generally two or three of these requirements play a 
pivotal role in sizing the wing, engine and take-off weight. 
Instantaneous turn rate is directly related to the amount of 
lift the wing produces. Sustained turn rate depends on lift 
and the amount of associated drag that engine thrust must 
overcome. The specific engine and wing sizes chosen to 
fulfil these maneuverability criteria will then affect the 
overall weight of the aircraft needed to meet mission radius 
requirements. Generally, wing and engine sizes must be as 
small as possible in order to minimize take-off weight. 

This makes it crucial for the designer to have sound 
aerodynamic data early in the conceptual design process to 
optimize performance and maneuverability for the smallest 
possible weight. The fact that the lift and drag of this 
aircraft are linked directly to the propulsion system 
magnifies the significance of being able to predict the 
aerodynamic characteristics with accuracy and convenience. 

As stated by the designer in Appendix A, page 63, the 
high aspect ratio nozzle can reduce the weight of the 
aircraft by improving the efficiency of the drag polar to a 
point where a wing with lower aspect ratio and/or area could 
be used. A smaller, lower aspect ratio wing generally 
weighs less, which thereby lowers the empty weight of the 
aircraft. The major question remaining is whether or not 





this weight savings can offset the weight increase of the 
installed nozzle. 

3 The Computational Model 
3.1 Pane 11ing 

The model as specified to LATJET was essentially a 
planar configuration with a "butterfly tail." Four major 
components, the forebody, canard, wing and v-tail modeled 



Note: V-TaH Omitted tor Ctartty 


Figure 6 b.lOJF-B Panel Model 


the aircraft, and a single component represented the jet 
sheet. (See Figure 6) All tips and sides have been inset 
1/4 of a panel width. The angled edge of the forebody has 
been approximated with a side parallel to the freestream. 
This was done in order to prevent extreme panel distortion 
caused by tapering a component to zero over a very small 
displacement in the Y direction. In order to represent the 
singularity at the jet, a 1/4 panel overlap was used as 
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recommended by reference 19. The jet sheet extends out 
three chord lengths past the trailing edge of the wing. 
Fredette supplied the NASA sharp edged transonic airfoil 
used in this study and seen in Figure 7. Since sizing the 
canard surface was one of the design objectives, the 
influence of wing camber was especially important in 


NASA TRANSONIC AIRFOIL W/SHARP L.E. 

t/c=.055 



x/c 


Figure 7 NASA Transonic, Sharp Leading Edged Airfoil 


determining the overall pitching moment. A camber line 
constructed with upper and lower surface coordinates was 
input to CAMBER.FOR. No twist, or camber variation was 
assumed in the spanwise direction. 
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3.2 Blowing Coefficient 

In calculating the blowing coefficients for this study, 
the nijVj term was assumed to be equal to the gross engine 
thrust at the specified flight condition. The mass flow 
rate and velocity across the nozzle exit plane were also 
assumed to be constant. This implies that the nozzle has a 
constant height side to side. This also means that 
varies in the spanwise direction proportionally to the 
chord. As far as flight parameters go, C„ was a function of 
Mach, altitude, and power setting. 

All engine data supplied by the designer was generated 
with a parametric engine deck employing the characteristics 
in Table II. The resulting gross thrust data was multiplied 
by an engine scale factor of 1.169. When calculating C M , 
the mjVj terms used were the gross thrust divided by the six 
panels across the exit plane of the nozzle. 

Table II b.lOJF-B Engine Characteristics 

Bypass Ratio 1.0 

Air Flow 350 lbm/sec 

Overall Pressure Ratio 30 

Turbine Inlet Temperature 3000°R 

To comprehend the effects of fanning the exhaust out, a 
reference case was created to see the characteristics of a 
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configuration with a smaller 2D vectoring nozzle. For the 
comparison case, the exhaust spanned over a width of 36 
inches (two panels) compared to 110 inches for the baseline. 

3.3 Nozzle Vector Angle 

The nozzle vector angle, or flap angle 6j was allowed 
to vary from -10 to +40 deg. It was assumed to act as a 
singularly blown flap, and have a constant mass flow rate in 
the spanwise direction. 

3.4 Modification for Control Surface Sizing 

To remove the influence of the existing canard for the 
canard sizing process, the canard component was terminated 
at the buttline even with the forebody. Individual moment 
curves were generated for this modified configuration. 
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IV. Verification 


1 Verification of Vortex Lat t ice _Code 

1.1 Two Dimensional Case With Camber (NACA 64-206) 

The verification studies were initiated by creating a 
flat plate model with a constant chord of one unit and a 
span of 100 units. This provided a pseudo 2D approximation. 
The lift curve slopes came out to 27t/rad, and a La( fO deg. 

Next, the functioning of the camber option was 
verified. Using the CAMBER program, a local angle of attack 
table was produced for a NACA 64-206 airfoil. This airfoil 
was chosen because it was relatively thin, its camber line 
could be easily generated by CAMBER.FOR, and its aerodynamic 
characteristics were readily available. The result are 
compared to values from Reference 1, in Table III. 

Table III LATJET Comparison for NACA 64-206 




LATJET 

Theory of Wing Sections 

C La 

(deg' 1 ) 

. 1095 

. 105 

a 1*0 

(deg) 

-2.07 

-1.5 

C Mac 


-.05 

- . 04 


It should be reemphasized that although the example airfoil 
has a 6% thickness, no thickness was modeled. 

1.2 Jet Sheet Convergence 

Convergence of the aerodynamic parameters with respect 
to jet sheet length was alluded to earlier. Jet sheets of 
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varying lengths were attached to the trailing edge of the 2D 
flat plate. C L was monitored as the sheet length increased. 
After approximately 3 chord lengths, the parameters converge 
to within 0.1% of the stabilized values. 

1.3 Three Dimensional Wing With Camber 

The b.lOJF-B wing was analyzed without its jet sheet. 

C La results were compared to analytic values. According to 
Hough (Ref 10), Eq (34) is accurate to 0.75% for 16>AR>1. 

C =. 2nAR 

^ AR+2.903+- 3 7-Z (34> 

AR 

The results are found in Table IV. Values for the wing with 
camber are also included. The negative zero lift angle of 
attack indicates that the camber is generating lift at a=0 
deg. 

Table IV Analytic and Computational Values for b.lOJF-B 
_ Wing _ 

Analytic Plain Wing Wing w/Camber 

C La (deg -1 ) .0431 .0463 .0463 

a L , 0 (deg) - 0 -.9365 


1.4 NADC Jet Flapped Fighter Study 

With the conventional aspects of LATJET checked out. 
the jet flap results needed evaluation. There is a great 
deal written about the jet flap. However, finding valid. 
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well presented wind tunnel data is very difficult. Quite a 
bit of data is available on blown flaps and upper surface 
blowing, but singularly blown jet flaps is another matter. 
Reference 9 was the only document located which involved a 
relatively simple configuration with a singularly blown 
flap. 

In reference 9, Hemmerly presents data for a family of 
tapered wings with and without tip plates and washout. The 
jet flap was intended strictly as a fixed high lift device 
which would vary its effect by modifying Cj. Air exits from 
a thin, spanwise, constant height slot placed just above a 
trailing edge with a small radius of curvature. A jet 
deflection angle of 65 deg was determined by tufting the 
trailing edge. Hemmerly made no mention as to how varying 
flight conditions and blowing coefficients affected the 
deflect ion. 

While reviewing data from this study, the importance of 
properly presenting all assumptions and book keeping schemes 
became apparent. In general, drag is assumed to be positive 
when measured in the direction of V x . With the reaction 
force of a jet flap to consider, the net force in this 
direction may actually be negative when measured by a 
balance in a wind tunnel. This would indicate that the 
thrust component of the jet has overcome the drag of the 
wing. This isn’t always obvious if the drag never becomes 
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negative, which was the case with the NADC study. The drag 
presented in the report included the thrust component of the 
jet flap (strictly as measured from the force balance). 

Drag curves compared here include the thrust of the jet. 

Another confusing aspect was the calculation of the 
blowing coefficient. The C u values as specified actually 
corresponded to a gross blowing coefficient, or Cj, as 
presented in Eq (9). Because the slot height was constant 
over the span, there was essentially a constant mass flow 
rate out the slot. All other things being equal, this meant 
that Cy varied over the span, inversely proportional to the 
chord. 

These facts were taken into account when a comparison 
was run using LATJET. The cases involve a simple tapered 
wing evaluated at three Cj values and an assumed flap 
setting of 65 deg. Results of the lift data can be seen in 
Figure 8. Lift coe fficient values in the three cases agree 
well with respect to both slopes and intercepts. 

Figure 9 shows the variation in drag coefficients and 
highlights the typical dilemma the vortex lattice method has 
with predicting drag. Each curve exhibits a parabolic 
relationship between lift and drag, but all three LATJET 
cases are optimistic compared to the wind tunnel results. 

It should be pointed out though, that the wind tunnel data 
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Figure 8 DTRC Fighter Comparison Lift Curves 


includes viscous effects whereas the computational data does 
not . 


LATJET Comparison W/DTRC Fighter 

S=12.04 f|\ cbar-1.834 ft, b-6.84 ft 



Figure 9 DTRC Fighter Comparison Drag Data 


Moment data in Figure 10 also shows similar trends. 
LATJET data appears to be slightly optimistic by under 
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predicting the magnitude of the pitching moment. However 
the vertical distribution of DTRC data for Cj=0 indicates 
that the reference position that the moment is calculated 
about may actually be the aerodynamic center, and not the 
c/4 as stated by Hemmerly. 


LATJET Comparison W/DTRC Fighter 

S=1 2.04 ft‘, cbar-1.834 ft, b=6.84 ft 
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Figure 10 DTRC Fighter Comparison Moment Data 
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V._Res ult s 

All cases were run from input files based on blowing 
coefficients. In the case of the configuration study, this 
refers to specific power settings. Each input file included 
individual cases to account for every combination of a and 

d j. 

1_F or c e Account i ng 

Before discussing the results themselves, the actual 
accounting of forces due to both the wing and jet should be 
explained. The lift coefficient. C L , as presented here 
includes both the lift of the wing and the component of jet 
thrust in the lift direction. In order to present a true 
powered drag polar, lift effects were not separated. Jet 
contributions to the lift coefficient can be calculated with 
the following equation. 

= C,sin(a+6,) (35) 

Since the drag is based on forces acting on each of the 
individual panels, the thrust of the jet must not be 
accounted for twice. It is popular convention for the 
designer to account for engine, and/or jet thrust 
independent from drag. This fact detaches the jet component 
from the system of aerodynamic forces acting on the 
aircraft. The induced drag presented here includes a near 




field calculation and a component due to the jet. C Djet . <' Djet 
is the result of the jet sheet never really reaching a 
horizontal position because of some finite value of downwash 
at (See Figure 11) This difference between theoretical 
and realized thrusts is accounted for as a component of 
drag. In terms of the LATJET output parameters. 



V 

c 

m 


Figure 11 C Djet Schemat ic 

= Cpi + CDJat + (36) 

With the sharp edged airfoil used on the b.lOJF-B, the 
Cpi value used included no leading edge suction. This 
equates to the situation of a sharp leading edge flap 
deflected directly into the relative wind. At high blowing 
coefficients, with increased circulation, the relative wind 
will have a definite vertical component with respect to the 


free stream. 







Table V lists the reference parameters used in 
calculating the nondimensiona1 coefficients shown for the 


b. 10JF-B. 

Table V b.lOJF-B Constant Parameters 


s ref (ft 2 ) 

481.5 

c (ft) 

19.36 

b ( ft ) 

37.58 

x cg (ft) 

31.6 


2_b .10JF-B Results 

Initial calculations showed that C„ values would fall 
within 2% of each other for identical power settings at both 
10,000 and 20.000 feet of altitude and a Mach number of 0.6. 
This indicates that for a given power setting, the change in 
engine gross thrust varies directly with density. Because 
of the similarity of the coefficients, only the data for 
20.000 feet is presented on the curves. The blowing 
conditions. BC, describe the current power setting and flap 
arrangement. Jet coefficient values corresponding to each 
blowing condition appear in Table VI. 

2.1 Influence of Blowing Coefficient 

Figures 12-17 display edited examples of data found in 
Appendix B, page 63. They primarily exhibit the effects of 
blowing coefficient on the various aerodynamic parameters 
plotted. Blowing conditions 3 and 4 have been removed from 
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Table VI b.lOJF-B Blowing Conditions 


BC 

Powe r Setting 

Cj 

Nozzle Width (in) 

0 

Unpowered 

0 

1 10 

1 

: Mil Powe r 

. 1658 

1 10 

2 

Max AB i 

. 2873 

110 

3 

Mil Powe r ; 

. 1658 

36 

4 

1 Max AB | 

. 2873 

! 36 


the plots to increase clarity. 

Figure 13 shows how the blowing coefficient increases 
the overall lift produced at a given angle of attack. At 
first, there appears to be no variation in the lift curve 
slope. Upon closer examination. Figure 12 reveals that C Lct 
actually varies linearly with respect to Cj. It increases 
approximately 10% over the range of Cj. 


Wright Lab b. 10JF-G 

M-0.6, c\-8°. 0% Leading Edge Suction 


0.057 - 



5[—10 
5j=0 


5j--20 


Sj-3C 


51=40 


Figure 12 


Variation in C la with Respect to C 










The variation of C L8j with respect to Cj also has a 
direct influence on lift. Figure 14 shows a non-linear 
relationship with respect to Cj. Since the study limited 


Wrign? I nn b. lO^F-B 

M=.6, 5 : =40°. 0'.’ I eudrg Suction 



Figure 13 Example b.lOJF-B Lift Data, 6j=40° 


itself to investigating blowing coefficients corresponding 
to only three power settings, the exact relationship is 
unc1 ear. 

The effect of Cj on the pitching moment is similar to 
its effect on the lift curve. As Cj increases, it shifts 
the curves downward. Figure 15. The slope of this curve has 
an important relationship to the static stability of an 
aircraft. The negative of its value is the static margin of 
the aircraft. 


SM = - 


Cm 

c L 


(37) 







Wright Lab b. 10JF-B 



0.001 

^ *CM)5 0.’ O.TT 0.2 ols ”^3 

CJ 

Figure 14 Variation of C L6j with Respect to Cj 

Wright Lab b.lOJF-B 

M=.6, 5i=4C 0 , O' 7 Lecding Edge Suction 
0.3 i-r-- tt--1-! 

! ■ 1 Lirocwer^b 



Figure 15 Example b.lOJF-B Moment Data, 6-40° 

Static margin is a convenient term which describes the 
distance of the CG from the aerodynamic center as a 
percentage of c. A negative static margin implies that the 
aircraft is statically unstable. Figure 15 shows that the 
static margin for the b.lOJF-B lies between -.21 and -.24 as 
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Cj varies from 0 to .2873. The relationship between SM and 
Cj is non- 1 inear. 
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Figure 16 Variation in Static Margin with Respect to Cj 


A mechanical flap analogy may explain the increase in 
static margin as Cj increases. Theoretically, a larger 
blowing coefficient represents a flap with more area. Since 
this increase in "area" is downstream of the trailing edge, 
the aerodynamic center shifts aft. This rearward shift in 
corresponds to the increasing static margin. 

The influence of Cj on induced drag can be seen in 
Figure 17. When compared to the drag polar for Cj=0, BC=0. 
the blown cases show induced drag penalties at C L =0. As «. 
or C L , increases, the blown curves begin to cross over the 
unpowered condition and show a decrease in induced drag. 

BC=1 crosses over at approximately C L =0.5 and BC=2 does so 
at about C L =.6. This demonstrates that an optimum Cj exists 
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for a particular flap deflection and lift condition. The 
induced drag increase at low a is due to the upset of the 
spanwise lift distribution from the flap deflection. 
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Figure 1~ Example b.lOJF-B Drag Data. 6^40“ 

2.2 Influence of Flap Deflection 

Example effects of flap deflection are presented in 
Figures 18-22. Each plot illustrates a blowing condition 
representing maximum afterburning thrust. BC=2. For 
clarity, flap deflections have been reduced to 0, 20. and 40 
degrees. 

Effects of 6j were found to be very similar to those of 
a mechanical flap. The lift curve in Figure 18, shows a 
definite lift increment for each level of 8_, with no 
discernable variation in C La . Looking at C La values with 
respect to 6j in Figure 19, a smoothly varying trend becomes 


apparent. 
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Figure 18 Example b.lOJF-B C L vs bj, BC=2 
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Figure 19 C La Variation with Respect to Oj 


The curves can be approximated with a cosine 
relationship. 

The rationale for a cosine approximation was that the 
projected "area" of the jet flap varies with respect to the 
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C ^ = . 00292 COSibj) * .051463 [deg' 1 ] for BC 1 

C^ = .00502 COS(6^ + .051463 [ deg l ] for BC 2 

cosine of the flap deflection. In Figure 19 the 
approximations overlap the actual values almost exactly. 

Induced drag in Figure 20 is minimized at 6j=0 and «=0. 
but there is a cross over phenomenon similar to the one seen 
in Figure 17. Knowledge of this is very beneficial since an 
aircraft will generally be flying at a fixed power setting 
and weight which would correspond to Cj and C L respectively. 
The optimized flap setting appear to reduce the induced drag 
by approximately 10%. 

The drag reduction may be short lived. This reduction 
comes with the corresponding increase in nose down pitching 
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Figure 20 Example b.lOJF-B Drag Data, BC=2 










moment, Figure 21. The control surface drag required to 
trim the moment to zero may be prohibitive. Again, like the 
lift data, there is a cosine variation in SM with respect to 
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Figure 21 Example b.lOJF-B Moment Data. BC=2 
Oj in Figure 22. The approximations are 
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Figure 22 Variation in Static Margin with Respect to 
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SM = . 01384 COSibj) - .24063 [deg' 1 ] 
SM = .02250 COS(6j) - .24063 [deg' 1 ] 


for BC 1 


for BC 2 

These values do not match quite as well at high flap 

2.3 Influence of Nozzle Width 

Figure 23 shows a lift curve with varying nozzle width. 
Both cases have the same value for Cj. 0.28"3. However, the 
narrow nozzle has a much larger value for C p across its 
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Figure 23 Example b.lOJF-B Nozzle Width Comparison, Lift 

center panels. With equal mass flow rates, the narrow 
nozzle must have a larger value of mj. There is no 
variation in slope. The increase in lift of the wider 
nozzle is the result of it acting on a much larger area of 
the wing. 








Figure 24 indicates a small induced drag increase due 
to the narrow nozzle. The overall shape of the polar is the 
same, and the is no cross over tendency. 
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Figure 24 Example b.lOJF-B Nozzle Width Comparison, Drag 


Pitching moments created by the narrow nozzle. Figure 25 
are less severe. 
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Figure 25 Example b.lOJF-B Nozzle Width Comparison, Moment 
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One assumption which needs to be addressed is jet 
thickness. As the jet width was reduced to 36 inches, the 
thickness increased to maintain a constant exit area. This 
case is presented as an example only. Any conclusions drawn 
from it must be regarded skeptically. Reference 12 presents 
a method for dealing with thick jets. 

All curves in Appendix B, except for the drag polars, 
show linear trends due to the absence of viscous effects. 

The CL vs a curves shows linear data up to and including 
ct=20 deg. This does not necessarily mean that this is an 
achievable, unstalled angle of attack, especially in the 
linear portion of the curve. C Lmax and a CL-ax values 
determined empirically will have to be faired in to the 
curves. 

3 Stability and Control Results 

Removal of the destabilizing canard changed the static 
margin to a positive value of approximately .02, Figure 26. 
This indicates aft movement of the aerodynamic center. The 
blown cases still show strong nose down moments which will 
have to be trimmed somehow. Moment variation due to Cj and 
5j is similar to the full configuration. Given the pitching 
moment without the canard, the size of the actual co.itrol 
surface can be calculated so that it produces an adequate 
level of control. Appendix C, page 86, contains the full 





set of moment graphs pertaining to the stability and control 
investigation. 
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Figure 26 Example b.lOJF-B With Canard Removed, Moment 

4 Planform Study Results 

Wing planform data can be found in Appendix D, page 92. 
Cases were run for each planform with variation of Cj and 6j 
about nominal values of 0.2 and 15 deg respectively. With 
respect to Sj, and Cj, the generic wing planforms behaved 
like the aircraft configuration. Similar trends were seen 
in all graphs. 

Since the primary purpose of this part of the study was 
to provide additional data to the designer on a variation in 
wing planform, the influence of AR and X characterized in 
the following example figures. 
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Figure 27 Lift Curve Variation with Respect to AR 


As the planform parameters, AR and X, varied, results 
typical of conventional wings occurred. Increasing AR 
raised the slope of the lift curve, Figure 27, and relaxed 
the drag polar, Figure 28. 
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Figure 28 Drag Polar Variation with Respect 
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The moment curve. Figure 29, shifted downwards and 
became steeper. This steeper, negative slope would have 
indicated an increase in the static margin had the moment 
been taken about a known or predicted CG instead of the 
planform quarter chord. 

Variation in \ showed little influence over the curves. 
There was essentially no change in the lift curve, 
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Figure 29 Moment Variation with Respect to AR 

Figure 30. The drag polar. Figure 31, also shows 
insignificant changes. There was however, an appreciable 
change in the pitching moment. Figure 32 shows the moment 
shifting downwards with the increase in A.. 
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L. ATJE T Wing Pianform Test Data 

AR=2.Q, CJ=0.2, DJ=15.G0, 0^ Leading Edge Suction 
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Figure 30 Variation of Lift with Respect to X 
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Figure 31 Variation in Drag with Respect to X 
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L AT JET Wing Planform Test Data 

AR-2.0. CJ-0.2, DJ-15.0C, 0% Leading Edge Suction 



Figure 32 Moment Variation with Respect to X 
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VI 


Conclusions and Recommendations 


The results of this study have shown unequivocal 
influence of the high aspect ratio nozzles on the 
aerodynamic characteristics of the Wright Lab b.lGJF-B 
fighter configuration. 

Raising the blowing coefficient showed positive shifts 
in the lift curve and increases in C La . Pitching moments 
moved in the nose down direction and the static margin 
increased. The drag polars showed that induced drag 
increases at low C L due to the disruption of the smooth 
spanwise lift distribution. As C L increased, the drag 
penalty disappeared, and the curves crossed over the unblown 
case, exhibiting an induced drag reduction. 

Variation of both the lift and moment curves with 
respect to 8j was very similar to the Cj cases. However, 
the slopes of both curves varied with respect to the cosine 
of 5j. The drag polar also behaved similarly with regarding 
8j. An induced drag penalty at low C L transitions to a 
benefit at high C L as the curve for a specific 6j crosses 
over the undeflected case. This suggests that an optimal 
flap setting exists for a given blowing condition and lift 
condition. These conditions correspond to the power setting 
and weight of the aircraft respectively. 
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Removal of the canard for the stability and control 
study resulted in a larger, positive static margin, and 
indicated a rearward shift in the aerodynamic center. 

All in all, the impact on the design can be summed up 
in the following points. The nozzle shows induced drag 
reduction at fixed power setting and optimum flap 
deflection. Drag reduction will be minimal at the partial 
power settings typical to cruise (Low Cj) and the relatively 
low C L values required. Trim drag as the result of 
countering the extreme pitching moment may overwhelm any 
drag reduction realized. 

The true benefits of the jet flap lie in the high lift 
region. This makes it most suitable for reducing take-off 
and landing distance and increasing sustained and 
instantaneous turn performance. However, this all hinges 
upon developing an effective way to counter the extreme 
pitching moments. 

Other aspects not included in this study must be 
accounted for prior to making any final decisions on design 
feasibility. Potentially, the high aspect ratio nozzle can 
reduce separation over the aft end. It may also cut down on 
the infra red signature of the exhaust plume. Nevertheless, 
production cost, weight and/or complexity may be 
prohibitive. The only complete answer will be the result of 
of a full blown multi-variate optimization study. 
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Appendix A: Ph ilos o phy of Designer 


The Wright Laboratory Configuration b.lOJF-B is a 
baseline layout for a future multi-role fighter design 
possibility. This single place, single engine concept 
incorporates a high aspect ratio jet flap type nozzle at the 
base of a wide, flat body. The design philosophy employed 
here is that the jet flap may upgrade the trimmed drag 
polars of low aspect ratio wings to the extent that its 
additional propulsion and integration weights can be offset 
by the reduced structural weight of the wings, canard and 
tails for a net reduction in weight empty, which will 
thereby improve the fuel load. The viability of this 
approach resides in the amount of aerodynamic improvement 
that can be derived from jet flap effects. 

Inspection of the 2 view drawing of the b.lOJF-B 
layout. Figure 1, and data sheet shows a design combining 
very high thrust and moderate wing loading. Highly 
favorable effects may make it possible to reduce both engine 
size and wing area in the interest of increasing fuel 
fraction. If the benefits allow changes in maneuvering wing 
loading to the region of 80 psf and a thrust-to-weight ratio 
of approximately 1.5, then the fuel fraction of a 30,000 lb 
fighter could be boosted 35% above the baseline case 
typified by the b.lOJF-B. 
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Table VII b.lOJF-B Dimensions and Weight 


Takeoff Gross Weight!lbs) 30,000 

Thrust to Weight Ratio 1.7066 

Wing Loading, takeoff (psf) 62.35 

Wing Loading, maneuver (psf) 57.00 

Overall Length (ft) 60.625 

Span (f t ) 31.573 

Maximum Body Width (ft) 14.543 

at FS 654.281 

Body Width at c/2 wing (ft) 13.413 

at FS 557.292 

Body Depth at FS 375 (ft) 5.333 

Overall Height, gear up (ft) 10.167 

Nominal C.G. at FS 451 


Raymond Fredette 
Senior Design Engineer 
WL/XPAD 


62 













Wright Lab b.lOJF-B 

M=.6, 5;=— 1 0°, 0% reading Edge Suction 



*(aeg) 


Figure 33 Lift Curve for 5j=-10° 
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Figure 34 Drag Polar for 6j=-10* 
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Wright Lab b. 10JF-B 
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Figure 35 Pitching Moment About the CG for 5j=-10° 
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Figure 36 Lift Variation Due to Cj for 6j=-l0° 
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Figure 37 Lift Curve for 6j=0° 
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Figure 38 Drag Polar for 8j=0° 
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Figure 41 Lift Curve for 6j=10° 
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Figure 44 Lift Variation Due to Cj for 6j=10° 
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Wright Lab b. 1 0 JF-B 
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Figure 45 Lift Curve for 8j=20° 
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Figure 47 Pitching Moment About the CG for 6j=20° 
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Figure 48 Lift Variation Due to Cj for 6j=20° 
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Figure 51 

Pitching Moment 

About the CG for 6j=30 
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Figure 52 Lift Variation Due to Cj for 6j=J0 
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Figure 53 Lift Curve for 5j=40° 
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Figure 54 Drag Polar for 6j=40 
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Figure 55 Pitching Moment About the CG for 6j=40° 
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Figure 56 Lift Variation Due to Cj for 6^40“ 
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Figure 57 Lift Curve for BC=0 
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Figure 58 Drag Polar for BC=0 
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Wright Lab b. 1 OJF-B 
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Figure 59 Pitching Moment About the CG for BC=0 
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Figure 60 Lift Variation Due to 6j for BC=0 
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Figure 63 Pitching Moment About the CG for BC=1 
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Figure 64 Lift Variation Due to 6j for BC=1 
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Figure 65 Lift Curve for BC=2 
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Figure 66 Drag Pola. for BC=2 
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CL 

Figure 67 Pitching Moment About the CG for BC=2 


Wright Lab b. 10 JF-B 

M=.6, BC=2, 0% Leading Edge Suction 



5i(degi 

Figure 68 Lift Variation Due to for BC=2 
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Wright Lab b. 10JF-B 

M=.6, 3C-3, 0 - Leading Edge Suction 
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Figure ^ 1 Pitching Moment About the CG for BC=3 


Wright Lab b.lOJF-B 

M=.6, BC=3, 0^ _eaa : ng Edge Suction 
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Figure 72 Lift Variation Due to Oj for BC=3 
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Wright Lab b. 1 OJF-B 

M=.6, 8C=4, O'" Leading Edge Suction 
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Figure 75 Pitching Moment About the CG for BC=4 


Wright Lab b. 1 OJF-B 

M=.6, BC=4. 0^ Leading Edge SucHon 
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Figure 76 Lift Variation Due to for BC=4 
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Wright Lab b. 10JF-B, NO CANARD, M=.5 

e>|=-10°. O ' Leading Edge Suct'on 

0.100 7 



Figure 77 Pitching Moment About CG, 6j=-10", No Canard 
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Figure 7g pitching Moment About CG, 6j-0°, No Canard 
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Wright Lab b. 10 JF-B. NO CANARD, M- 

5j=10°, 0^ Leading cage Suction 
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Figure ?9 Pitching Moment About CG, dj = 10°, i\o Canard 


Wright Lab b. 10JF B, NO CANARD, M=.6 

5j=20°, O 07 Leaa : ng Edge Suction 


Figure 80 Pitching Moment About CG, 6j=20°, No Canard 
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Wright Lab b. 10JF-B, NO CANARD, M=.5 

5|-30°, 0 U ' Leadrg Edge Suc K on 

0.C50 i - 1 -r—-r- 1 - —1 


x; - £ 


30 D i-r-r- i ■ ■ -r-P >"-“i i T-r- , t i . ) — 

-0.50 -0.00 0.5C 1.0C 1.50 


Figure 81 Pitching Moment About CG, 6j=30®, No Canard 


Wright Lab b. 10JF-B, NO CANARD, M=.6 

5j=40°, 0% Leading Edge Suction 
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Figure 82 Pitching Moment About CG, 6j=40°, No Canard 
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Wright Lab b.lOJF-B, NO CANARD. M-.6 

CJ^O.OOOC, 0“' Leading Edge Suction 
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Figure 83 Pitching Moment About CG, BC=0, No Canard 


Wright Lab b. 10JF-B, NO CANARD, M=.6 

CJ=0.1658, 0% Leading Edge Suction 
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Figure 84 Pitching Moment About CG, BC=1, No Canard 
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Wright Lab MOJF-B, NO CANARD, M=.5 

CJ-0.28 7 3, O'- Leading idge Sweden 
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Figure 85 Pitching Moment About CG, BC= 


No Canard 











LATJET Wing Planfcrm Test Data 

AP-2.0, CJ=0.2, T R'-0.05, 0 ' J Leading Edge Suction 



• 0.5 *-1 r i- : 4 
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c<(oeg,» 


Figure 86 Lift Curve Comparison for 6j, AR=2, A. = 0.0 


LaTJET Wing Planform Test Data 

AR=2.0, DJ-'S 0. iP—0.05. 0 ” Leading Edge Suction 
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Figure 87 Lift Curve Comparison for BC. AR=2. A=0.0 











L.ATJET Wing Pcnform Test Data 

AR=2.0, CJ=0.2, TR-0.05, 0" Leading Eoge Suc‘ : on 
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Figure 88 Drag Polar Comparison for 0 j, AR=2, X=0.05 


LATJET Wing Planform Test Data 

AF-2.0, DJ=15.0, TR-0.05, 0" Leading Edge S^ctio^ 
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Figure 89 Drag Polar Comparison for BC, 


AR= 2 . 


A = 0.05 


94 








CMc/4 


LATJET Wing Planform Test Data 

AR=2.0, CJ=0.2, TR=0.05, Cr 5 ' Leading Ecge Sucrcn 
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Figure 90 Moment Comparison for 6j, AR=2. X=0.05 


LATJET Wing Planform Test Data 

AR=2.0, DJ=15.0, TR=0.05, 0% Leading Edge Suction 
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Figure 91 Moment Comparison for BC. AR=2, X=0.05 
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LATJET Wing Planform Test Date 

AR-2.0, CJ=0.2, TR-0.05, O ’ 7 Leaa’ng Edge Suction 
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Figure 92 Lift vs 6 jf AR=2. A. = 0.05 


LATJET Wing Planform Test Data 

AR=2.0, DJ=15.0, T R-0.05, 0% Leading Edge Suc+ion 
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Figure 93 Lift vs Cj, AR=2, A. = 0.05 
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LATJET Wing Planfcrm Test Data 

AP-2.0, CJ-0.2, T P-0.25. O l " r Leading Edge Section 
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Figure 94 Lift Curve Comparison for Oj, AR=2. X=0.25 


LATJET Wing D lanform Test Data 

AR=2.0, DJ=15.0, TR=0.25, 0^ Leading Edge Suc+ion 
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Figure 95 Lift Curve Comparison for BC, AR=2, X=0.25 
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L AT JET Wing Picnforrm Test Data 

AR=2.0, CJ-0.2, 7R-0.25, O ’ Leading Ecge $uc T 'on 

2 1- 1 -T-—T- t I 
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Figure 96 Drag Poiar Comparison for Oj, AR=2. >1 = 0.2 


LATJET Wing Planform Test Data 

AR=2.0, DJ=15.0, TR-0.25. 0d Lecding Edge Sjcticr 
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Figure 97 Drag Polar Comparison for BC, AR=2, X = 0.2 
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LATJET Wing Planform Test Data 

AF=2.0, CJ=0.2, TR=0.25, C l >r Leading Edge Suci'on 
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Figure 98 Moment Comparison for 6j. AR = 2, X=0.25 


LATJET Wing Planform Test Data 

AR-2.0, DJ=15.0, TR-0.25, 0% Leading Edge Sucncr 
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Figure 99 Moment Comparison for BC. AR = 2, A. = 0.25 
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LATJET Wing D lanform Test Date 

AR-2.0, CJ=0.2, T3-0.25, 0 v Leading Edge Suction 
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Figure 100 Lift vs Oj, AR=2, X=0.25 


LATJET Wing Planform Test Data 

AR=2.0, Dj= 15.0, TR=0.25, 0^ Leading Edge Suctior 
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Figure 101 Lift vs Cj, AR=2, X=0.25 
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LATJET Wing Planfcrm Test Data 

AP=3.0, CJ=0.2, T P-0.05, 0% Leading Edge Suc^on 
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Figure 102 Lift Curve Comparison for Oj, AR = 3. A, = 0.05 


LATJET Wing Planform Test Data 

AP=3.0, DJ= 1 5.0, TR=0.Q5. 0*3 Leading Edge Suction 
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Figure 103 Lift Curve Comparison for BC, AR = 3, >1 = 0.05 
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L AT JET Wing P’cnform Test Data 

AR=3.0, CJ=0.2, TR-3.05, CT Leading Edge Suc^on 
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^0.00 0.05 0.10 0.15 0.20 0.25 0.50 0.55 0.40 

Figure 104 Drag Polar Comparison for Oj, AR = 3. >1 = 0.05 


LATJET Wing Planform est Data 

AR-3.0, DJ-15.C, TR=0.05, Leading Edge Suction 
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Figure 105 Drag Polar Comparison for BC, AR=3, A=0.05 
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LATJET Wing Planfo^m Test Data 

AP=5.0, CJ=0.2, TR-0.05, 0 ” .eccinc Eoge Sucron 
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Figure 106 Moment Comparison for oj, AR=3. X=0.05 


LATJET Wing Plantorm Test Data 

AR=3.0, DJ=15.0, TR-0.G5, 0 Jr _eaoing Edge Su-ricr 
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Figure 107 Moment Comparison for BC, AR = 3, A. = 0.05 
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LATJET Wing °lanform Test Date 


AR-i.O, CJ-0.2, TR-0.C5. 0 " _eaa'ng Eage Suclio’ 
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Figure 108 Lift vs Oj, AR = 3, X=0.05 


LATJET Wing Planform Test Data 

AR-3.G, DJ-15.0, tr-c. 05, 0^ Leading Edge Suc’icr 
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Figure 109 Lift vs Cj. AR=3. *=0.05 











LATJET Wing Planfcrm Test Data 

AP=3.0, C J-0.2, T( ?-0.25, O'" Leading Edge Sudan 
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Figure 110 Lift Curve Comparison for 6j, AR=3. A = 0. 


LATjET Wing Planform Test Data 

AP=3.0, DJ=15.0, TR-0.25, 0- Leading Edge Section 





O.OC 2.GO 4.0C 6.0C 8.00 "Q.OG 12.0C 14.00 16.00 

o(aeg) 
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LATJET Wing Pcnform Test Data 

AR-3.0, CJ-0.2. TR-C.25. Q -' Leading Edge $u*: f 'cn 



gure 112 Drag Polar Comparison for Oj, AR = 3. A. = 0 . 


LATJET Wing Planform Test Data 

AR-3.C, DJ— 15.0, TR^C.25, 0^ ; _eoding Edge Suction 


CJ-0 



gure 113 Drag Polar Comparison for BC. AR=3, X=0. 
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L AT JET Wing Planfo^m Test Data 

AP=3.0, CJ=0.2, TR=0.25, 0" .eoding Edge Sucron 
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Figure 114 foment Comparison for 5j, AR=3. A. = 0.25 


LATJE^ Wing Plantorm Test Data 

AF-3.0, DJ-15.0, TR=C.25, 3^ Leading Edge Suction 
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Figure 115 Moment Comparison for BC. AR=3. X=0.25 
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LATJET Wing Planform Test Date 

AR-3.0, CJ—0.2, TR-0.25, 0'3 Lead'ng Edge Sucilo n 
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Figure 116 Lift vs 6j, AR=3, X = 0. 25 


LATJET Wing Plantorm Test Data 

AR=3.0, DJ-15.0, TR-C.25. 0% Leading Edge Suction 



Figure 11^ Lift vs Cj. AR=3, X=0.25 
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l_ATJET Wing Planfcrm Test Data 

AP-3.0, CJ=0.2, T R-0.50, 0"’ Lecaing Edge Suc + ion 
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Figure 118 Lift Curve Comparison for 8j, AR=3, X=0.5 


LATJET Wing Planform Test Data 

AR=3.0, DJ=15.0, TR=0.50, 0^ Leading Edge Suction 
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Figure 119 Lift Curve Comparison for BC, AR=3. X=0.5 
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LATJET Wing PiGnform Test Data 

AR-3.0, CJ=0.2, 7R=0.50, Q% Leading Edge Suc t; on 
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Figure 120 Drag Polar Comparison for 5 r AR=3, X=0.5 


LATJET Wing Planform T est Data 

AR-3.0, DJ-15.0, TR=C.50, 0^ Leading Edge Suction 
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Figure 121 Drag Polar Comparison for BC, AR=3, X=0.5 
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L AT JET Wing Planform Test Data 

AP=3.0, CJ=0.2, TR-0.5C, 0' r _eacing Ecge Sucron 
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Figure 122 Moment Comparison for oj. AR=3, A. = 0.5 


LATJE7 Wing Plantorm Test Data 

AR-3.0, DJ-15.0, TR=0.5C, 0” Leading Edge $uc*icr 
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Figure 123 Moment Comparison for BC. AR=3. A. = 0.5 








LATJET Wing D !anform Test Date 

AR-3.0, CJ-0.2, IR-C.50, 0" Leaa : ng Zcge SuCi,'o~ 
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Figure 124 Lift vs Oj, AR = 3, X = 0.5 


LATJET Wing Planform Test Data 
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Figure 125 Lift vs Cj, AR=3 
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LATJET Wing Planfcrm Test Data 

AR-5.0, CJ=0.2, T R-0.50, 0 > Leading Edge Suction 
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Figure 126 Lift Curve Comparison for 6j, AR=5. i=0.5 


LATjET Wing Planform Test Data 

AR=5.0, 0J= 1 5.0, tR=O.50, 0% Leading Edge Suction 
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Figure 127 Lift Curve Comparison for BC, AR=5, X=0.5 
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LATJET Wing PiGnform Test Data 

AR=5.0, CJ=0.2, 7R=0.50, Leading Edge Suc^on 
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Figure 128 Drag Polar Comparison for 8j, AR=5, X=0.5 


LATJET Wing Planform Test Data 

AR=5.0, DJ=15.0, TR=0.50, 0 % Leading Edge Suction 
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Figure 129 Drag Polar Comparison for BC, AR=5, X=0.5 
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LATJET Wing Planform Test Data 

AP=5.0, CJ=0.2, TP=0.50, 0 v Leading Edge Sucron 
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Figure 130 Moment Comparison for 6j, AR=5, A,=0.5 


LATJET Wing Planform Test Data 

AF=5.0, DJ=15.0, TR^O.50, 0% Leading Edge Suction 
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Figure 131 Moment Comparison for BC, AR=5, X=0.5 
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LATJET Wing Planform Test Date 

AR=5.0, CJ=0.2, TR-0.50, 0- Leaa : ng Ecge Suction 
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Figure 132 Lift vs 6j, AR=5, >1=0.5 


LATJET Wing Planform Test Data 

AR-5.0, DJ-15.0, tr-q. 50, 0% Leading Edge Suction 



Figure 133 Lift vs Cj, AR=5. A.=0.5 
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LATJET Wing Pianfcrm Test Data 

AP-5.0, CJ=0.2, TC ?-1.00, 0 ur Leading Edge Sue ■ ion 
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F'gure 134 Lift Curve 

Comparison 

for 6j 

AR=5 , k= 1 


LATJET Wing Planform Test Data 

AF=5.0, DJ=15.0, TR= 1 .0C, 0*2 Leading Edge Suction 
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Figure 135 Lift Curve Comparison for BC. AR = 5. >1=1 
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Figure 136 Drag Polar Comparison for AR=5, X=1 


LATJET Wing Flanforrm Test Data 

AF=5.0, DJ=15.0, TR».t..0C. On Leading Ecge Suction 



0.00 0.05 0.10 0.15 0.20 C.25 0.30 0.35 0.40 0.45 0.50 
CL 


Figure 137 Drag Polar Comparison for BC, AR=5. X=1 
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LATJET Wing Planform Test Data 
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Figure 138 Moment Comparison for 6j, AR=5, X=1 
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Figure 139 Moment Comparison for BC, AR=5, A.= 1 
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Figure 140 Lift vs 6 j, AR=5, X=1 
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Figure 141 Lift vs C J( AR=5 , X=1 
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